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ABSTRACT 

The  subject  of  this  report  is  the  mechanical  behavior  of  stiffened  plates,  basic  structural 
components  of  ships  and  submarines.  The  buckling  loads  and  modes  of  grillages  subjected 
to  axial  compression  without  and  with  lateral  pressure  are  calculated  using  finite  element 
based  eigenvalue  analyses. 
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INTRODU  CTION 


Stiffened  plates  are  basic  structural  components  of  ships  and  submarines.  In  our  earlier 
report*,  we  have  calculated  the  axial  buckling  loads  of  grillages  with  the  use  of  a  well-known 
finite  element  code.  In  the  present  work,  we  continue  the  calculation  of  buckling  loads  for 

new  grillages. 

The  grillages  now  modeled  consist  of  3  base  plates,  2  longitudinal  edge  plates,  4  longi¬ 
tudinal  stiffeners,  and  2  transverse  stiffeners.  Two  sets  of  grillages  (I  and  II)  having  different 
plating  thicknesses  and  stiffener  dimensions  are  considered  (Figures  1,  2  and  12,  13).  We 
define  the  following  grillage  dimensions,  of  which  the  first  4  are  common  to  all  grillages; 

al  =  length  between  transverse  stiffeners  =  96'' 

a2  =  length  between  transverse  stiffener  and  grillage  end  =  96" 

bl  =  width  between  longitudinal  stiffeners  =  27" 

b2  =  width  between  outer  longitudinal  stiffener  and  grillage  side  =  7.5" 
b3  =  width  of  longitudinal  edge  plates 

tl  =  thickness  of  inner  base  plate  and  longitudinal  edge  plates 

t2  =  thickness  of  outer  base  plates 

dwl  =  depth  of  longitudinal  webs 

twl  =  thickness  of  longitudinal  webs 

dfl  =  depth  of  longitudinal  flanges 

tfl  =  thickness  of  longitudinal  flanges 

(iw2  =  depth  of  transverse  webs 

tw2  =  thickness  of  transverse  webs 

df2  =  depth  of  transverse  flanges 

tf2  =  thickness  of  transverse  flanges 

*D.  A.  Danielson  and  D.  P.  Kihl,  “Buckling  of  Ship  Grillages,”  Naval  Postgraduate  School  Technical 
Report  NPS-MA-96-002,  September  1996. 
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In  these  grillages,  the  outer  base  plates  also  have  thickness  tl  in  regions  of  width 
adjacent  to  the  transverse  stiffeners.  The  material  is  isotropic  steel  with  Young’s  modulus 
£;  =  3  X  10’^  psi  and  Poisson’s  ratio  v  =  .3.  For  grillages  I,  the  base  plates  and  longitudinal 
edge  plates  are  HSLA  steel  with  yield  stress  85  ksi,  whereas  the  longitudinal  and  transverse 
stiffeners  are  HS  steel  with  yield  stress  60  ksi.  For  grillages  II,  the  steel  types  are  interchanged; 
the  plates  are  HS  steel,  whereas  the  stiffeners  are  HSLA  steel. 

All  grillages  have  the  same  imposed  boundary  conditions: 

(1)  One  end  of  a  grillage  has  all  3  displacement  components  zero  and  all  3  rotation 
components  zero. 

(2)  The  other  end  of  the  grillage  where  the  force  is  applied  has  uniform  axial  displacement 

with  the  other  2  displacement  components  zero  and  all  3  rotation  components  zero. 

(3)  The  ends  of  the  transverse  stiffeners  have  vertical  displacement  zero  but  the  other  2 
displacement  components  and  all  3  rotation  components  are  free. 

(4)  For  the  pressure  cases,  the  longitudinal  edge  plates  are  also  restrained  vertically  at 
10.5"  intervals  along  the  outer  sides. 

(5)  All  other  nodes  on  the  plates  edges  are  completely  free. 

To  each  of  these  grillages  we  apply  compressive  axial  loading  together  with  one  of  3 
normal  loads;  no  pressure,  an  initial  pressure  loading  of  10  psi  (simulating  fluid  on  the 
plating  side),  or  an  initial  pressure  loading  of  -10  psi  (simulating  fluid  on  the  stiffener  side). 
Furthermore,  to  assess  the  effects  of  the  longitudinal  edge  plates,  we  consider  both  uncut 
and  cut  edge  plates.  This  makes  a  total  of  12  different  cases  to  study. 

The  finite  element  code  used  is  MSC  Nastran/Patran.  The  pre-processing  and  post¬ 
processing  is  done  with  Patran  (version  6.2).  The  buckling  analysis  (solution  105)  is  done 
with  Nastran  (version  69).  The  meshing  is  done  entirely  with  Quad  4  plate  elements.  Note 
that  this  is  a  full  finite  element  model  of  the  entire  grillage  and  does  not  assume  symme¬ 
try.  We  use  linear  geometry  and  linear  material  properties  for  all  cases.  Computations  are 
performed  on  a  Silicon  Graphics  Indy  Workstation. 
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GRILLAGES  I* 

The  first  set  of  grillages  has  the  following  dimensions: 
b3  =  6" 
tl  =  .4375" 
t2  =  .5" 
dwl  =  6" 
twl  =  .17" 
dfl  =  4" 
tfl  =  .215" 
dw2  =  14" 
tw2  =  .255" 
df2  =  5" 
tf2  =  .42" 

The  mesh  size  used  is: 

mesh  length  of  all  elements,  mesh  width  of  base  plate  elements  1.5 
mesh  width  of  web  and  edge  plate  elements  =  1" 
mesh  width  of  flange  elements  =  .5" 


♦Subsequent  to  the  completion  of  the  grillages  I  cases,  it  was  discovered  that  the  longitudinal  edge  plates 
should  have  extended  4"  below  the  base  plating  (instead  of  4"  on  the  stiffeners’  side).  For  the  unpressurized 
uncut  case,  the  buckling  force  is  then  1860  kips  (as  opposed  to  1865  kips).  Since  the  difference  between  the 
two  results  is  so  small,  we  did  not  correct  the  grillages  I  cases.  (In  the  grillages  II  cases,  the  longitudinal  edge 
plates  are  exactly  centered  along  the  base  plating.)  In  addition,  in  the  positively  pressurized  grillages  I  cases, 
the  longitudinal  edge  plates  should  have  been  constrained  vertically  along  the  base  plating’s  side  rather  than 
the  stiffeners’  side.  For  the  positively  pressurized  cut  case,  the  buckling  force  is  then  1350  kips  (  as  opposed 
to  1352  kips).  Since  the  difference  between  the  two  results  is  so  small,  and  to  ensure  the  end  displacement 
under  pressure  -10  psi  is  the  negative  of  the  end  displacement  under  pressure  10  psi,  we  did  not  change  the 
side  along  which  the  constraints  are  applied.  (In  the  pressurized  grillages  II  cases,  the  longitudinal  edge 
plates  are  constrained  vertically  along  the  base  plating’s  side.) 
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NO  PRESSURE 


First  we  consider  the  linear  response  under  an  axial  compressive  force  of  1865  kips.  The 
prebuckling  state  is  one  of  uniform  axial  compression  with  little  bending.  The  magnitude 
of  the  maximum  stress  occurs  in  the  thinner  regions  of  width  6"  in  the  outer  base  plates 
(Figure  3).  Note  that  the  magnitude  of  the  maximum  stress  is  about  50  ksi,  considerably 
less  than  the  yield  stress. 

The  bifurcation  buckling  force  is  1865  kips.  The  buckling  mode  is  a  deformation  in¬ 
volving  primarily  bending  with  little  stretching  (Figure  4).  The  inner  base  plate  buckles  into 
a  square  quilt  of  half-wavelength  24".  The  adjacent  webs  also  buckle  with  half-wavelength 
24".  The  flanges  mainly  just  rotate  about  their  center  line.  The  edge  plates  also  just  rotate 
about  their  line  of  attachment  to  the  base  plate.  In  our  earlier  report  we  called  this  buckling 
mode  TRIPPING. 

To  assess  the  effects  of  the  longitudinal  edge  plates,  we  next  cut  the  element  sides  lying 
along  vertical  lines  at  intervals  10.5".  The  force  required  to  buckle  this  grillage  with  48  edge 
cuts  is  1721  kips  (8%  less  than  the  uncut  case).  The  buckling  mode  is  similar  but  with 
relatively  larger  deformations  near  the  cut  edges  (Figure  5). 

POSITIVE  PRESSURE  (Fluid  on  Plating  Side) 

Next  we  consider  the  linear  response  under  a  normal  pressure  loading  of  10  psi  and  an 
axial  compressive  force  of  1436  kips.  The  prebuckling  state  now  involves  axial  compression 
and  bowing  out  of  the  entire  grillage  with  lateral  torsional  bending  of  the  stiffeners  (Figure 
6).  However,  the  magnitude  of  the  stress  is  still  less  than  the  yield  stress,  except  at  the  ends 
of  the  central  longitudinal  stiffeners. 

In  the  presence  of  positive  pressure,  the  bifurcation  buckling  force  is  reduced  to  1436  kips 
(23%  less  than  the  unpressurized  case).  The  buckling  mode  is  similar  to  the  unpressurized 
case  (Figure  7). 

To  again  assess  the  effects  of  the  longitudinal  edges,  we  again  make  48  vertical  cuts  in 
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the  edge  plates.  This  cut  grillage  buckles  at  an  axial  compressive  force  of  1352  kips  (6%  less 
than  the  uncut  pressurized  case).  Since  the  edge  plates  are  now  restrained  vertically  along 
their  upper  edges  at  each  of  the  cuts,  the  buckling  mode  does  not  show  as  much  deformation 
near  the  cut  edges  as  the  cut  unpressurized  case  (Figure  8). 

NEGATIVE  PRESSURE  (Fluid  on  stiffener  Side) 

Finally  we  consider  the  linear  response  under  a  normal  pressure  loading  of  -10  psi  and 
an  axial  compressive  force  of  1516  kips.  Note  that  the  prebuckling  compressive  stress  in  the 
center  of  the  longitudinal  stiffeners  is  now  larger  than  in  the  base  plates  (Figure  9).  The 
magnitude  of  the  stress  is  now  less  than  the  yield  stress  everywhere. 

In  the  presence  of  negative  pressure,  the  bifurcation  buckling  force  is  1516  kips  (19% 
less  than  the  unpressurized  case).  The  buckling  mode  now  has  half- wavelength  48  (Figure 

10). 

The  cut  grillage  under  negative  pressure  buckles  at  an  axial  compressive  force  of  1427 
kips  (Figure  11)  (again  6%  less  than  the  uncut  pressurized  case). 
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ship4:  Buckling  force1865  kips  -  Magnitude  of  displacement 
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ship12:  Pressure  10  psi,  Force  1436  kips  -  Von  Mises  stress 
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Pressure  10  psi,  Buckling  force  1436  kips  -  Magnitude  of  displacement 
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ship17:  Pressure  -10  psi,  Force  1516  kips  -  Von  Mises  stress 


ship16:  Pressure  -10  psi,  Buckling  force  1516  kips  -  Magnitude  of  displacement 
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8:  Pressure  -10  psi,  Buckling  force  1427  kips  -  Magnitude  of  displacemen 


GRILLAGES  II 

The  second  set  of  grillages  has  the  following  dimensions: 
b3  =  8.5" 
tl  =  .375" 
t2  =  .4375" 
dwl  =  5.5" 
twl  =  .3125" 
dfl  =  4.25" 
tfl  =  .3125" 
dw2  =  13.5" 
tw2  =  .25" 
df2  =  5" 
tf2  =  .4375" 

The  mesh  size  now  used  is: 

mesh  length  of  all  elements,  mesh  width  of  base  plate  elements  =  1.5" 

mesh  width  of  edge  plate  elements  =  .85  . 

mesh  width  of  webs  =  5". 

mesh  width  of  longitudinal  flanges  =  .53125 

mesh  width  of  transverse  flanges  =  .625" 
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NO  PRESSURE 


The  magnitude  of  the  maximum  stress  in  the  prebuckling  state  under  an  axial  force  of 
1552  kips  is  about  35  ksi,  much  less  than  the  yield  stress. 

The  bifurcation  buckling  force  is  1552  kips  for  the  uncut  grillage  and  1406  kips  for  the 
cut  grillage.  The  buckling  modes  are  similar  to  the  corresponding  grillages  I  cases  (Figures 
14  and  15). 

POSITIVE  PRESSURE  (Fluid  on  Plating  Side) 

The  magnitude  of  the  maximum  stress  in  the  prebuckling  state  under  a  normal  pressure 
of  10  psi  and  an  axial  compressive  force  of  1403  kips  is  about  46  ksi,  still  considerably  less 
than  the  yield  stress. 

In  the  presence  of  positive  pressure,  the  bifurcation  buckling  force  is  1403  kips  for  the 
uncut  grillage  and  1293  kips  for  the  cut  grillage.  The  buckling  modes  are  similar  to  the 
corresponding  grillages  I  cases  (Figures  16  and  17). 

NEGATIVE  PRESSURE  (Fluid  on  stiffener  Side) 

In  the  presence  of  negative  pressure,  the  bifurcation  buckling  force  is  1685  kips  for  the 
uncut  grillage  and  1555  kips  for  the  cut  grillage.  The  buckling  mode  now  has  half-wavelength 
of  about  20"  (Figures  18  and  19). 
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ship21:  Buckling  force  1552  kips  -  Magnitude  of  displacement 
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Ship23:  Pressure  10  psi,  Buckling  force  1403  kips 
Magnitude  of  displacement 
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ship24:  Pressure  10  psi,  Buckling  force  1293  kips 
}  Magnitude  of  displacement 


mmmmm 


fMWSS’S 

Wi 

V- 

m 


m 


m 

i/SSmm 


i 


////W/ 


‘^/iir 

wir»|| 

I 


Af 


JW 


H 

iisiitr 


mSm^^Wmrnmi 

mP^ifussssi 


25 


CONCLUSIONS 


GRILLAGES  I 


Continuous  Edges 

Cut  Edges 

Pressure  0 

-.3047",  1865  kips 

-.2963"  ,  1721  kips 

Force  0 

,001643",  10  psi 

.001758",  10  psi 

TABLE  1:  END  DISPLACEMENT,  FORCE  OR  PRESSURE 


Continuous  Edges 

Cut  Edges 

Pressure  0 

1865  kips  ,  1 

1721  kips  ,  .923 

Pressure  10  psi 

1436  kips  ,  .770 

1352  kips  ,  .725 

Pressure  -10  psi 

1516  kips  ,  .813 

1427  kips  ,  .765 

TABLE  2;  COMPRESSIVE  BUCKING  FORCE,  NORMALIZED  BUCKLING  FORCE 
GRILLAGES  II 


Continuous  Edges 

Cut  Edges 

Pressinre  0 

-.2542",  1552  kips 

-.2477"  ,  1406  kips 

Force  0 

.001938"  ,  10  psi 

.002130"  ,  10  psi 

TABLE  3:  END  DISPLACEMENT,  FORCE  OR  PRESSURE 


Continuous  Edges 

Cut  Edges 

Pressure  0 

1552  kips  ,  1 

1406  kips  ,  .906 

Pressure  10  psi 

1403  kips  ,  .904 

1293  kips  ,  .833 

Pressure  -10  psi 

1685  kips  ,  1.086 

1555  kips  ,  1.002 

TABLE  4:  COMPRESSIVE  BUCKING  FORCE,  NORMALIZED  BUCKLING  FORCE 
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The  axial  displacement  at  the  end  of  a  grillage  under  the  compressive  buckling  force 
or  the  positive  pressure  is  given  in  Tables  1  or  3.  Since  the  prebuckling  state  is  linear,  the 
axial  displacement  under  both  loads  is  the  linear  superposition  of  the  two  displacements  in 
a  column.  To  get  the  displacement  at  any  other  load,  we  simply  multiply  by  the  displace¬ 
ment/load  ratio  obtained  from  Tables  1  or  3. 

The  predicted  buckling  loads  for  the  12  cases  are  summarized  in  Tables  2  and  4.  Dividing 
each  force  by  the  upper  left  entry  in  the  table,  we  obtain  the  normalized  buckling  force,  from 
which  we  can  immediately  see  the  reduction  in  buckling  force  due  to  applying  pressure  or/ and 
cutting  the  edge  plates. 

Judging  from  our  earlier  work,  the  buckling  loads  we  have  now  calculated  are  likely  to  be 
greater  than  collapse  loads  to  be  measmed  in  future  experiments.  Unmodeled  factors  which 
would  tend  to  decrease  the  buckling  loads  include  geometric  imperfections,  residual  stress, 
and  relaxed  boundary  conditions.  On  the  other  hand,  we  have  calculated  only  the  forces 
required  to  initiate  buckling,  but  the  grillages  are  likely  to  have  significant  postbuckling 
strength  before  ultimate  collapse. 
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